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Abstract 

The Thorium abundance for the red giant COS 82 in the Ursa Minor dwarf spheroidal galaxy is 
determined based on a high resolution spectrum. This is the first detection of actinides in an extra Galactic 
object. A detailed abundance pattern is determined for 12 other neutron-capture elements from the atomic 
number 39 to 68. These elements are significantly over-abundant with respect to other metals like Fc (> 
1 dex) and their abundance pattern agrees well with those of the r-process-enhanced, very metal-poor stars 
known in the Galactic halo, while the metallicity of this object ([Fe/H]~ —1.5) is much higher than these 
field stars ([Fe/H] ~ —3.0). The results indicate that the mechanism and the astrophysical site that are 
responsible for neutron-capture elements in COS 82 is similar to that for field r-process-enhanced stars, 
while the condition of low mass star formation is quite different. An estimate of the age of this object 
based on the Th abundance ratio is discussed. 

Key words: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: individual 
(Ursa Minor COS 82) — galaxies: dwarf 
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1. Introduction 

The rapid process (r-process) is known to be responsi- 
ble for about half of the abundances of elements heavier 
than Fc in solar-system material. Since the reaction oc- 
curs through very neutron-rich, unstable nuclei, the un- 
derstanding of the process is still a challenging issue in 
nuclear physics. The astrophysical sites of this process 
are still unclear, though numerical simulations have been 
made for many possibilities (e.g. Wanajo & Ishimaru 
2006). 

An important progress in this field in the past decade 
is the discoveries of very metal-poor stars having large 
over-abundances of r-process elements. Such objects are 
expected to record the yields of nucleosynthetic events in 
the early Galaxy including the r-process. Several objects 
having very large excesses of r-process elements (e.g. Eu) 
with respect to other metals ([Eu/Fc]> 1) arc known in 
the Galactic halo (Sneden et al. 2003; Hill et al. 2002; 
Christlicb et al. 2004; Barklem et al. 2005), and are called 
as 'r-IF stars (Beers & Christlieb 2005). A remarkable 
result for these stars is that the abundance patterns of 
neutron-capture elements are very similar. Objects having 
such large excesses of r-process elements are only found at 
the metallicity around [Fc/H]= —3.0, and no object having 
[Eu/Fe]> +1.0 is known in the range of [Fc/H]> —2.5 in 
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the Galactic halo, though the reason for such a metallicity 
distribution is still unknown. 

However, the situation is quite different in some dwarf 
spheroidal galaxies (dSph), in which high abundance ra- 
tios of r-process elements are found even in stars with 
relatively high metallicity. The chemical abundance stud- 
ies for individual stars in the local group dSph made large 
progresses in the past several years, thanks to high reso- 
lution spectrographs mounted at 8-10 m telescopes (e.g. 
Shetrone et al. 2001, 2003; Venn et al. 2004). One 
surprising result derived by these abundance studies is 
the existence of stars showing extremely large enhance- 
ments of heavy neutron-capture elements with relatively 
high metallicity. The most remarkable star is the red gi- 
ant COS 82 in the Ursa Minor dSph. High resolution 
spectroscopy of this object was first obtained by Shetrone 
et al. (2001; in their paper, this star is referred to as 
'199'), who reported that this star is moderately metal- 
poor ([Fc/H]^ —1.5) but shows large enhancement of 
heavy neutron-capture elements ([Eu/Fe]=:1.49). Our pre- 
vious observing obtained a red spectrum of this star with 
higher quality, and has successfully detected some other 
heavy elements (e.g. Dy, Er), confirming that the abun- 
dance pattern of elements heavier than Ba almost com- 
pletely agrees with the r-process abundance pattern esti- 
mated from solar-system material (Sadakane et al. 2004). 
This result is surprising, because no moderately metal- 
poor star ([Fe/H]> -2) in our Galaxy is known to date to 
have such large excesses of r-process elements. Tsujimoto 
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& Shigeyama (2002) applied their supernova-driven star 
formation model and proposed that this phenomenon can 
be explained by assuming a large velocity dispersion of 
the gas from which low mass stars like COS 82 formed. 

In this Letter, we report the detection of the Th ab- 
sorption line at 5989 A for this object by the re-analysis 
of our high resolution spectrum obtained by Sadakane et 
al. (2004). The actinide Th has an isotope whose half 
life is 14.05 Gyr, and is produced only by the r-process. 
The production of actinides is very sensitive to the envi- 
ronment of the sites (e.g. entropy, timescale of the ex- 
plosion) according to model calculations for the r-process. 
Therefore, the Th abundance can be a very strong con- 
straint on the study of the origin of the neutron-capture 
elements in this object. While Th abundance has been 
determined for more than ten field metal-poor stars and 
globular cluster stars, our measurement for COS 82 is the 
first determination of an actinide abundance for an extra 
Galacitic object. We also discuss the possibility to apply 
the Th abundance ratio to estimate the age of this object. 

2. Observation, Analysis and Results 

A high resolution spectrum of COS 82 was obtained by 
Sadakane ct al. (2004) with the Subaru Telescope High 
Dispersion Spectrograph (HDS, Noguchi et al. 2002). The 
spectrum covers 4400-7100 A with a resolving power of 
40,000. The signal-to-noise ratio per 1.8 km -1 pixel is 60 
at 5900 A. Standard data reduction procedures are carried 
out with the IRAF cchclle package 1 as described by Aoki 
et al. (2005). Equivalent widths for isolated absorption 
lines are measured by fitting gaussian profiles. 

A standard abundance analysis was made for the mea- 
sured equivalent widths using the model atmospheres of 
Kurucz (1993). The analysis was made primarily based 
on the line list of Sadakane et al. (2004). We added 
some weak Fe I lines, which enable one to determine the 
micro-turbulence more accurately. The effective temper- 
ature (T c ff) of 4300 K, determined by Sadakane et al. 
(2004), was adopted in the analysis. The micro-turbulence 
(^turb) and gravity (g) are determined so that the derived 
Fe abundance is not dependent on the strengths of Fe I 
lines nor on the ionization stage, respectively. The re- 
sults (logg = 0.6 dex and fturb = 1-7 km s _1 ) agree with 
those of Sadakane et al. (2004) within the estimated er- 
rors. The metallicity is assumed to be the Fe abundance 
([Fc/H]). The Fe abundance derived by the present anal- 
ysis is [Fe/H] = —1.42, adopting the solar abundance of 
loge(Fe) = 7.45 (Asplund et al. 2005). The results of the 
abundance analysis are given in table 1. 

The Th abundance is determined from the Th II 5989 A 
line. The oscillator strengths of Th lines, including the 
5989 A line, are obtained by Nilsson et al. (2002). This 
line was investigated by Yushchcnko et al. (2005) for the 
field metal-poor star HD 221170. The line was measured 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 



for the same object, using a higher quality spectrum, by 
Ivans et al. (2006), who reported that the abundance de- 
termined from this line agrees with those from other Th 
II lines. We performed an analysis for HD 221170 us- 
ing the spectrum obtained with the Subaru Telescope and 
the model atmosphere (Kurucz 1993) with the parame- 
ters adopted by Ivans et al. (2006), and confirmed that the 
derived Th abundance agrees very well with their result. 
We note that the partition function of the Th II was cal- 
culated using the energy levels listed by Blaise & Wyart 
(1992). 

Though telluric absorption lines, probably due to wa- 
ter vapor, exist in this wavelength range, we found that 
the Th II 5989 A line of COS 82 is not distinctively af- 
fected by them. Figure 1 shows the observed spectrum 
along with the synthetic ones for three Th abundances. 
The data for spectral lines other than the Th II line were 
adopted from the list of Kurucz & Bell (1995). The wave- 
length of the observed spectrum is determined from 137 
Fe I lines in the spectrum. We identified the absorption 
feature at 5989.3 A as a Nd II line. Although a Nd II 
line at 5989.378 A with the lower excitation potential 
of 6005.270 cm" 1 is listed by Kurucz & Bell (1995), J. 
E. Lawler (private communication) suggests that the line 
should be the transition from 19758.540 to 3066.755 enr 1 
at A air = 5989.312 A. Since the transition probability of 
this line is unknown, we assume it to be logg/ = —2.05. 
This g/-value and the Nd abundance determined from 
other spectral regions well reproduce the strength of the 
line for COS 82 as well as for HD 221170. 

The Th abundance is determined by x 2 -htting of syn- 
thetic spectra to the observed one from 5988.9 to 5989.2 A, 
which is not severely affected by the 5989.3 A feature. The 
derived abundance is loge(Th)= —0.25. The 2a level of 
the fitting error is 0.07 dex. The continuum level is esti- 
mated for the spectrum around the absorption features. 
The uncertainty of the continuum placement (~ 0.5%) re- 
sults in an possible error of 0.03 dex in the Th abundance. 
We estimate the line broadening for relatively clean lines 
in the spectrum, and found that the line widths are pri- 
marily determined by the instrumental broadening, and 
no significant broadening by macro-turbulence nor rota- 
tion is detected. The possible errors in the Th abundance 
due to the adopted broadening parameter (±0.5 km s -1 ) 
is 0.03 dex. 

The abundances of other elements are determined using 
the line list of Sadakane et al. (2004) with some modifi- 
cations. The transition probabilities of Nd II, Sm II and 
Gd II are updated, adding several new lines to the list, 
using the data of Den Hartog ct al. (2003), Lawler ct al. 
(2006) and Den Hartog et al. (2006), respectively. The 
effects of hyperfine splitting for Ba, La, Eu, and Pr are in- 
cluded using the line data obtained by McWilliam (1998), 
Lawler, Bonvallet, Sneden (2001), Lawler et al. (2001), 
and Ginibre (1989), respectively. We also measure the 
O abundance from the [O I] 6300.3 A line. The results 
are listed in table 1. An upper limit of the Os abun- 
dance is derived from the Os I 5584 A line. Estimates 
of upper limits for other neutron-capture elements are 
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also attempted, but no meaningful result is derived from 
our spectrum. The abundances of Fe and most neutron- 
capture elements derived by the present analysis agrees 
with the results of Sadakane et al. (2004) within the er- 
rors (see below). Exceptions are Eu and Gd abundances, 
for which the differences between the two measurements 
are about 0.25 dex. This would be due to the difference of 
the spectral features and line data adopted in the analysis. 
The abundances of Fe and some neutron-capture elements 
determined by Shetrone et al. (2001) agree well with our 
results, while their Ce and Sm abundances are more than 
0.3 dex higher than ours. We suspect that their results 
are uncertain because of the low S/N ratio in the blue 
region of their spectrum where Ce and Sm lines used in 
their analysis exist. 

The abundance uncertainties for elements for which 
more than 10 lines are measured are estimated by ctA" 1 / 2 , 
where a is the standard deviation of the abundances from 
individual lines and N is the number of lines used. A 
typical uncertainty of the equivalent width measurement 
from our spectrum is 6 mA. Wc performed abundance 
analyses for Fe I lines with equivalent widths from 30 to 
200 mA, changing the equivalent widths by 6 mA. The 
effects of the change of equivalent widths on the derived 
abundances are 0.15-0.20 dex, which agree with the a 
(0.18 dex) for Fc I (<tfci)- For elements whose abun- 
dances are determined from less than 10 lines, the ran- 
dom error is estimated by OFeiA -1 / 2 . Errors due to the 
uncertainty of atmospheric parameters are estimated by 
calculating abundances changing atmospheric parameters 
by A(T off )=150 K, A(log 5 )=0.3 dex, A([Fe/H])=0.2 dex, 
and A(uturb) = 0.3 km s _1 . These errors are added, in 
quadrature, to the random errors estimated above, and 
are given in table 1 (cr to tai)- We found, however, that 
the abundances of the neutron-capture elements deter- 
mined from the ionized species systematically scales by 
changes of the adopted atmospheric parameters. For the 
discussion based on the abundance pattern, we calculate 
the average of the abundance changes for heavy neutron- 
capture elements (La-Er) by changing the atmospheric 
parameters, and adopt the deviation from this systematic 
abundance change for individual elements as the errors 
due to the uncertainties of atmospheric parameters. This 
results in quite small errors (< 0.05 dex). These errors 
and the random errors are added in quadrature, and are 
also given in table 1 as cr n _ 

(see text). 
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Fig. 1. Comparisons of synthetic spectra for the region 
including the Th II 5989 A line with the observed spec- 
trum. The Th abundances assumed in the calculations are 
loge(Th)= -0.15,-0.25,-0.35 (solid lines), and -co (dotted 
line). 












3. Discussion and concluding remarks 
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The upper panel of figure 2 shows the abundance 
patterns of neutron-capture elements for COS 82, the 
solar-system r-process component (Simmerer et al. 2004), 
and the two r-process-enhanced field stars CS 22892- 
052 (Sncdcn ct al. 2003) and CS 31082-001 (Hill et al. 
2002). The abundances are normalized to the average 
of Ba Er ones. As reported by Sadakane et al. (2004), 
the overall abundance pattern of heavy neutron-capture 
elements (Ba-Er) in COS 82 agrees well with the r- 



Fig. 2. Upper panel: Scaled abundances of neutron-capture 
elements for COS 82 (filled circles), CS 22892-052 (open cir- 
cles), CS 31082-001 (crosses), and the solar-system r-process 
component (solid line). The abundances arc normalized to the 
average of the abundances of the 10 elements from Ba to Er 
in the logarithmic scale. Lower panel: Abundance differences 
in the logarithmic scale between COS 82 and CS 22892-052 
(open circles) and that between COS 82 and the solar-system 
r-process component (filled circles). Solid lines indicate the 
averages of the abundance differences for Ba-Er. 
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Table 1. Elemental Abundances Results 



species 


l0g£sun 


loge 


[X/Fej 


N 


"total 


^n — cap 


Fc I 


7.45 


6.03 




137 


0.27 




Fe II 


7.45 


6.02 




13 


0.18 




I 


8.70 


7.45 


0.17 


1 


0.2 




Y II 


2.21 


1.22 


0.42 


6 


0.23 


0.14 


Zr II 


2.59 


1.82 


0.65 


1 


0.24 


0.19 


Ba II 


2.17 


1.36 


0.60 


3 


0.30 


0.21 


La II 


1.13 


0.52 


0.81 


8 


0.16 


0.06 


Ce II 


1.63 


0.76 


0.55 


4 


0.20 


0.14 


Pr II 


0.80 


0.20 


0.82 


7 


0.19 


0.11 


IN a 11 


1 A K 
1.40 


1 nn 

1 .UU 


u.y { 


41 


U.l { 


U.UO 


Sm II 


0.98 


0.90 


1.34 


2 


0.21 


0.13 


Eu II 


0.52 


0.34 


1.24 


3 


0.17 


0.11 


Gd II 


1.09 


1.10 


1.43 


3 


0.16 


0.12 


Dy II 


1.17 


1.13 


1.38 


3 


0.17 


0.10 


Er II 


0.97 


0.73 


1.18 


3 


0.17 


0.11 


Th II 


0.09 


-0.25 


1.08 


1 


0.15 


0.08 



process abundance pattern. In the lower panel, the abun- 
dance differences between COS 82 and the solar-system 
r-process component, as well as that between COS 82 and 
CS 22892-052, are shown. The Th abundance of COS 82 
is lower than the normalized value of the solar-system r- 
proccss component, while that agrees very well with the 
value of CS 22892-052. This point is discussed below in 
more detail. The abundance of the light neutron-capture 
element Y in COS 82 is also lower than the normalized 
value in the solar-system r-process component as found in 
field r-process-enhanced stars. 

Although such small deviation exists, an important re- 
sult derived from the comparison is that the abundance 
pattern of neutron-capture elements, including Th, in 
COS 82 agrees with the pure r-process abundance pattern. 
According to model calculations for the r-process, the ac- 
tinidc abundances with respect to other neutron-capture 
elements are very sensitive to the model parameters (e.g. 
electron fraction, entropy; Wanajo et al. 2002). The 
agreement of the Th abundance ratio to other neutron- 
capture elements in COS 82 with that of CS 22892-052 
indicates that the mechanism and astrophysical sites that 
are responsible for the heavy elements in COS 82 is com- 
mon with those for the r-process-enhanced field stars. 
Recent abundance studies for r-process-enhanced stars 
suggest that the r-process occurs under a quite limited 
condition (e.g. supernova explosions whose progenitors 
have mass in a very narrow range). Our result shows that 
this is also applied to stars even in dSph. 

The metallicity of COS 82 is more than one order 
of magnitude higher than those in the field r-II stars 
([Fc/H]~ —3). The chemical composition of field stars 
with such metallicity is believed to be affected by a num- 
ber of nucleosynthesis events. However, the metallicity, as 
well as the amount of neutron-capture elements, should be 
related to how the yields from a nucleosynthesis event are 
captured by the clouds from which next generation, low 
mass stars formed, and the neutron-capture elements of 



COS 82 might be provided by a single event. For instance, 
Tsujimoto & Shigcyama (2002) proposed that the large 
excesses of neutron-capture elements in a few objects in 
dSph, including COS 82, are the result of a larger veloc- 
ity dispersion of the gas from which low mass stars like 
COS 82 formed in dSph than in the Galactic halo. That is, 
the yields of one supernova are recorded in the next gen- 
erations of stars with smaller mixing with metal-deficient 
(or metal-free) gas in dSph. 

Finally, we discuss an estimate of the age of COS 82 
from the abundance ratio between Th and other stable el- 
ements, which could be an independent calibration for the 
age measurements for this galaxy. If the age of this star 
is shown to be quite long (> 10 Gyr), the star formation 
of this dSph should be completed with a short time scale. 
Th has a half life of 14.05 Gyr, and the age of the object 
(more strictly, the duration from the r-process event) is 
estimated from the Th abundance ratio to other stable 
elements, if its initial abundance ratio is given. Here the 
average of abundances of the nine elements from La to Er 
is used as the reference 2 (E sta bie)- The Th/E sta bie ratio 
of COS 82 agrees very well with that of CS 22892-052, 
while that is 0.3 dex lower than the value in the r-process 
component in solar-system material (figure 2). 

Unfortunately, the abundance ratio of Th/E sta bic pro- 
duced by the r-process is not constant. Indeed, the ac- 
tinides including Th are significantly enhanced, with re- 
spect to other neutron-capture elements, in the r-II star 
CS 31082-001, compared to that of the solar-system r- 
process component (Hill et al. 2002). However, studies 
of Th abundances for a number of r-process enhanced 
field stars show that the abundance ratios distribute in 
a rather narrow range (—0.6 <Th/E sta bic < — 0.1), and 
CS 22892-052 is one of the objects that have the low- 

2 Ba is excluded from the reference elements, because its abun- 
dance error is relatively large due to the strengths of the lines 
used in the analysis. The conclusion in this paragraph is, how- 
ever, unchanged if this element is included in the calculation. 
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est ratio. Assuming that the initial abundance ratio of 
Th/E s tabic for COS 82 is within this range, and that the 
field halo stars including CS 22892-052 are quite old (their 
ages are around 12-13 Gyr), COS 82 is also an old object. 
(If a higher value of the initial Th/E sta bic ratio is assumed, 
the derived age becomes longer). 

The error (an— cap) m the determination of the Th 
abundance ratio is 0.08 dex, while the scatter of the 
abundances of the above nine stable elements around 
those of the solar-system r-process component, as well 
as of CS 22892-052, is 0.11 dex. We estimate the error 
of the Th/E s tabio to be 0.13 dex, which results in the 
uncertainty of the age estimate to be 6 Gyr. Hence, 
though the object is suggested to be as old as field halo 
stars, the age estimate using Th chronometer gives only a 
weak constraint. A measurement of U/Th ratio, which is 
a more sensitive chronometer, for this object will enable 
one to estimate the age with a much smaller error (< 2 
Gyr). 

We are grateful to Dr. J. E. Lawler for providing useful 
Nd II line data. 
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